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ABSTRACT 

The interaction of open and closed field lines at coronal hole boundaries is widely accepted to be 
due to interchange magnetic reconnection. To date, it is unclear how the boundaries vary on short 
timescales and at what velocity this occurs. Here, we describe an automated boundary tracking 
method used to determine coronal hole boundary displacements on short timescales. The bound- 
ary displacements were found to be isotropic and to have typical expansion/contraction speeds of 
<2 km s~^, which indicate magnetic reconnection rates of <3 x 10~^. The observed displacements 
were used in conjunction with the interchange reconnection model to derive typical diffusion coeffi- 
cients of <3 X 10^^ cm^ s~^. These results are consistent with an interchange reconnection process 
in the low corona driven by the random granular motion of open and closed fields in the photosphere. 
Subject headings: Sun: corona — Sun: evolution — Magnetic fields 
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1. INTRODUCTION 

Coronal holes (CHs) are regions of the solar atmo- 
sphere with reduced emission and density and hence 
appear as da rk areas at X-ra.y and extreme ultraviolet 
wavelengths (Ahschulcr et ajj 119721 : iVaiana et aD 119761 : 
Ic'h ap ma.n fc B romage ,2002i). As regions of predomi- 
nantly open magnetic field, CHs are known to be dom- 
inated by a single polarity. This characteristic is of- 
ten use d to distinguish CHs from other low-intensity 
regions (Harvev & Recelyl |2002 IScholl fc Habball 120081 : 
[Krista & Gallagher 200°^. Here, we study the small- 
scale changes in the CH boundaries that are thought to 
result from magnetic reconnection occurring between CH 
open fields and small clo sed fields found in the q uiet Sun 
(QS) as weh as in CHs (|Madiarska et al.ll2004[ ). We in- 
vestigate the driving force behind the open and closed 
field interaction and determine the magnetic reconnec- 
tion rate and the diffusion coefficient based on the ob- 
served displacement in CH boundaries. 

A number of authors have previously studied the 
large-scale and long-term evolution of CHs and es- 
tablished the changes in the total CH area, fiux, 
the latitudinal distribution over the solar cycle and 



terchange recon n ection" (IWang fc SheelevI 
Wang et al .' 1996; Wa ng fc SheelevI 120041: IFisia 



Sun (iBravo & Stewartl 


19971 iChaoman & Bromagell2002l: 


iHarvev & Recelvl 12002 


: IWang et all 119961: IWand 120091). 



Although the small-scale changes in CH boundaries have 
only been studied in d etail in recent years using high- 
resolution observations (iMadiarska fc Wie gelmannl l2009l : 
ISubramanian et al.ll2010t) . it was alreadv suggested in the 
late 1970s that the small-scale and short-term evolution 
of CH boun daries is strongly connected to bright point 
occurrences (jNolte et al.lll978l: lDavislll985l ) and that re- 
connection could be a means of flux-transportation and 
"diffusion" of the open flux (,Nolte et al. 1978). A more 
in-depth discussion on reco nnection occurring at CH 
boundaries was provided by iMuUan fc AhmadI ()1982D . 
who proposed that reconnection was also a driver of the 
solar wind emanating from CHs. 

The process of continuous reconnection at CH 
boundaries later became better known as "in- 
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Raiu et al.i ,2005: Edniondson et al.l 120101) . As argued 
bv iFisk fc SchwadronI (120011 ) and iFisM (l20Q5l) . the 
transfer of magnetic flux along the solar surface can be 
considered as a diffusion-convection process. Despite 
all of the processes involved being convective in nature, 
the non-stationary random supergranular convection is 
included in the diffusion term while those that are uni- 
form (such as differential rotation and meridional flow) 
are included in the convection term. As supergranular 
motions are too slow for flux transportation, interchange 
reconnection was suggested as a faster mechanism for 
field-line diffusion. Interchange reconnection allows for 
the transfer of magnetic fiux from one point to another 
through magnetic r econnection betwee n open and closed 
field lines, which iFisk fc SchwadronI (|2001[ ) modelled 
using a diffusion-convection equation, 



dt 



(1) 



where Bz is the radial magnetic field, u is the uniform 
convective flow velocity and 



(5rf 
2St 
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is the diffusion coefficient that depends solely on the open 
field displacement distance (Sr) and the time elapsed 
{5t). A graphical representation of interchange recon- 
nection is given in Figure 1, where a QS loop borders the 
CH open magnetic field. Magnetic reconnection between 
the open and closed field leads to the displacement of 
the open field and thus the CH boundary moves further 
out. During this process the larger loop is replaced by 
a smaller loop with oppositely oriented footpoint polari- 
ties imbedded inside the CH, while the open field line is 
"relocated" to the former foot-point of the larger loop. 

Interchange reconnection at CH boundaries is also 
supported by observational results, e.g., the preva- 
lence of very small loops inside CHs with heights 
<10 Mm and larger loops with heights 10-20 Mm 
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Figure 1. Interchange reconnection: The location of the CH and 
the QS are indicated on the cartoon with a Une representing the 
boundary between the two regions. The scenario on the left shows a 
QS loop located near a CH open field which yields magnetic recon- 
nection due to the proximity of the opposite polarity loop footpoint 
to the open field. The result of the reconnection is shown on the 
right - a small loop is created inside the CH with opposite mag- 
netic orientation and the open field line "relocated" from within 
the CH to the further foot-point of the previously larger loop. This 
process can ex plain the rig id rotation as well as the growth of CHs 
(modified from rFisH[2005f) . 

outside CHs HWiegelmann fc So lank? '2004"), the pres- 
ence of bidirec t ional jets (Madiarska ct al. 20Q3; 
IDovle et all 120061: IKamio et all 120091V and the bound- 
ary displacements observed due to the emergence 
and disappearance of bright small-scale loops in the 
form of bright points (iMadiarska fc WiegelmannI 120091 : 
iSubramanian et all I2010D . While the observed small- 
scale dynamic processes appear to be in agreement with 
the interchange recon nection model, it was argued by 
lAntiochos et al.l ()2007[ ) that intermixed open and closed 
fields of the interchange reconnection model would lead 
to the formation of currents sheets, which cannot be 
maintained in a low plasma-beta corona. Instead, the 
quasi-steady models (e.g. potential field source surface 
and magnetohydrodynamic models) explain the long- 
term evolution of the global open magnetic flux by the 
diffusion of active region magnetic flux a nd its trans- 
port to high latitudes by the meridional flow fWan g et al.l 
120001 12010). Also, the quasi-steady models suggest that 
open field regions of the same polarity are connected by 
"corridors" , since disconnected open fields (such as those 
created i n the interchange reconnection process) are no t 
allowed (lAntiochos et a"Lll2007HEdmondson et al.ll2009[ ). 
lAntiochos et al.l ()2007[) suggested that the "corridors" 
connecting the open field regions could be wide contain- 
ing highly dynamic, intermixed open and closed fields 
which could accommodate interchange reconnection. 

In this Letter, we study the short-term evolution of 
CH boundaries by objectively determining CH bound- 
aries based on intensity thresholding methods, as well 
as introducing a new automated boundary tracking tool 
that allows us to follow the changes that occur (detailed 
in Section|2]). In SectionjSjthe magnetic reconnection rate 
and diffusion coefficient are determined using the inter- 
change reconnection model, and the measured boundary 
displacement velocities. The processes that influence the 
CH boundary evolution are discussed in Section |4l 

2. OBSERVATIONS AND DATA ANALYSIS 

Finding objective CH boundary intensity thresholds 
is fundamental in deriving consistent physical prop- 
erties (e.g., area, magnetic field strength) throughout 
the solar cycle. Many authors have previously cho- 
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Figure 2. CH contours in an EIT 195 A Lambert projection map 
(a) and the corresponding sub-image (b). MDI sub-image with the 
CH contours (c) (>50 G fields: yellow, <-50 G fields: blue) and 
the corresponding -Blos histogram (d) (magenta: CH, black: QS). 
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sen arbitrar y intensity thresho l ds to define CH bound- 
aries ( e.g.. Ifiaker et alj [20071: [Abramenko et alj 120091: 



iMadiarska fc WiegelmannI 12009 ). I n order to improye 
the lo calisation of CH boundaries, iKrista fc Gallagheil 
([2009') developed an automated intensity thresholding al- 
gorithm to determine the boundary threshold between 
CHs and the QS based on local intensity histograms. 
The same method is used in this study for a smaller 
region of interest, a 31° x 35° box centered on each 
CH studied. Here we use 12 minute cadence 195 A im- 
ages from t he SOHO Extreme ultravio let Imaging Tele- 
scope fEIT: [Delaboudiniere et al.lll995[ ) to locate low in- 
tensity regions, and 96 minute cadence magnetograms 
from the Michelson Doppler Imager (MDI: [Scherrer et al.l 
119951 ) to detect predominant polarities and hence iden- 
tify CHs. The magnetograms were also used to determine 
the magnetic properties of the CHs and the surrounding 
QS. For this study, eight relatively small CHs were cho- 
sen arbitrarily between 2008 September 8 and October 
17, and tracked for two to six days. The CHs studied 
included both predominantly positive and negative po- 
larity holes, with equal numbers located in the north- 
ern and southern hemispheres. All of the CHs studied 
were non-fragmenting and of relatively simple morphol- 
ogy throughout their evolution. 

The KIT images were initially transformed into Lam- 
bert equal-area projection maps to aid the thresholding 
algorithm and to simplify the tracking of the region of in- 
terest. Figure [2^1. shows the KIT Lambert projection map 
with the CH contours and the box around one of the CHs 
studied (only the largest CH in the box was a candidate 
in our study). Figures [2)d and[2t show the KIT and MDI 
close-up images of the box with the CH contours, while 
Figure [2Ji shows the histogram of the line-of-sight mag- 
netic field (-Blos) from the largest CH (magenta line) and 
the surrounding QS (black line) excluding small satellite 
CHs. The skewness of the i3LOS histogram is used to 
establish whether a predominant polarity is found, thus 
allowing the differentiation of CHs (skewness > 0.5) from 
other low intensity regions (skewness < 0.5). The agree- 
ment between the visual and the analytically-obtained 
CH boundary, together with the predominant polarity 
found in the contoured regions, validate the appropriate- 
ness of the boundary positions and justifies the CH status 
of the detected regions. In principle, the intensity of a 
CH boundary can change due to the line-of-sight obscu- 
ration at higher longitudes and latitudes. However, all 
CH candidates in this study lie within ±50° in longitude 
and latitude, which lowers both the obscuration effect in 
EUV images and the errors in line-of-sight magnetic field 
measurements . 

Once a CH was identified, consecutive boundary ob- 
servations were compared. Considering an example of 
contour pairs as shown in Figure [3^, the centroid of the 
earlier observed contour (solid line, observed on 2008 
September 8 at 05:48 UT) was considered as the origin 
of the frame of reference for both contours. The later 
contour (dashed line, observed on 2008 September 8 at 
17:00 UT) will have moved with the synodic solar rota- 
tion rate compared to the earlier contour, due to solar 
rotation. For this reason a correction was applied to 
the second contour coordinates. The locations of both 
contours were measured radially from the earlier cen- 
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Figure 3. (a) The CH contour observed at 05:48 UT (solid line) 
and 17:00 UT (dashed line) on 2008 September 8. The frame of 
reference is the earlier contour centroid and angle is measured anti- 
clockwise from the right direction of the earlier contour centroid. 
(b) The radial location of the boundary plotted as a function of 
position angle (solid: earlier; dashed: later), (c) Boundary dis- 
placement velocities as a function of position angle, (d) Boundary 
displacement velocity distributions as a function of position angle 
for all observations of the CH. The dotted lines in (a), (b), and (c) 
mark the boundary locations at 0°, 90°, 180° and 270°. The error 
bars in (d) correspond to ±1 a for each 20° bin. 

troid in 5° segments and anti-clockwise from west of the 
earlier contour's centroid (Figure [5)d). Each boundary 
segment location was determined as the median of the 
contour point distances from the reference centroid. The 
displacement of a segment (Sr) and the time difference 
between the two observations yielded the displacement 
velocity of the segment. This calculation was repeated 
to get all the boundary displacement velocities over 360° 
(Figure [3j;) . The displacement velocities measured for 
the whole observation period of a single CH were also 
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determined, as shown in the 2D histogram of Figure 
Here, dark regions correspond to more frequent occur- 
rences of the boundary velocity values, where the ob- 
servations were binned in 20° and 0.4 km s~^ intervals. 
These calculations were further extended to include ev- 
ery boundary displacement velocity of each CH in the 
study, as detailed in the following Section. 

3. RESULTS 

3.1. Coronal hole boundary displacement velocities 

In the case of the aforementioned example CH, simi- 
lar average velocities were observed over all angle space 
(Figure [3jl) . This indicates that there is no preferred di- 
rection in the boundary evolution of the CH, and hence 
the boundary evolution is isotropic. The value of the av- 
erage signed velocities, (v), is ~0 km s~^ in all angle bins, 
which indicates that there is no consistent expansion or 
contraction in the CH during the time of observation. 
Furthermore, the majority of ±1 a values for all angle 
bins are below 2 km s^^ (as is the case for all eight CHs 
studied here) , which is in good agreement with the aver- 
age velocities of random granular motions. 

In order to extract characteristic velocities, we combine 
all CH observations over all angles and consider the un- 
signed boundary velocity distribution (Figure H^). From 
this distribution we determine typical unsigned velocities 
to be <2 km s~^. The positive and negative velocities 
tend to cancel when averaged over time, giving an aver- 
age signed velocity of ~0 km s~^. However, this does not 
imply that the CHs studied showed no variation in size 
during the observing period, rather that the expansion 
and contraction were on comparable scales. 

3.2. Magnetic reconnection rates 

Magnetic reconnection between CH open fields and 
small-scale QS loops are expected to occur in the tran- 
sition region and low corona. Here we can make the as- 
sumption that as a small-scale loop reconnects with the 
open magnetic field in the transition region it will conse- 
quently cause a displacement between open and closed 
field lines in the inner corona. In fact, the displace- 
ment at the transition region should be roughly pro- 
portional to the displacement detected in EUV images 
if the surrounding QS has a relatively stable magnetic 
configuration. The Alfven speed, wa, was calculated us- 
ing VA = B/ y/ATTnipne, where B is the magnetic field 
strength and n,, is the electron number density. Typ- 
ical values for the magnetic field {B — 2 G) and den- 
sity (rte = 10^ cm^'^) in the lower co rona were used 
from lSchriiver fc van BallegooiienI ()2005[ ) to calculate the 
Alfven speed, which was found to be ~600 km s~^. The 
magnetic reconnection rate was determined from. 
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(3) 



where Vin is the infiow speed, the relative speed of the 
opposite polarity magnetic field lines being brought to- 
gether prior to reconnection. The distribution of dis- 
placement velocities was used for v-m and substituted to- 
gether with the Alfven speed into Equation |3] to obtain 
the reconnection rate distribution. This distribution is 
shown in Figure |41d, where typical reconnection rates are 
<3 X 10-3. 



3.3. Magnetic diffusion coefficients 

We calculate the magnetic diffusion coefficient for each 
5° segment of each CH observation using Equation [2] 
by substituting the boundary displacement distances in 
the segment {6r) and the time elapsed (St) between 
consecutive observations. The resulting distribution of 
diffusion coefficients is shown in Figure |4j;, with val- 
ues typically <3 x 10^^ cm^ s'-^. Using typical loop 
footpoint distances and ch a racter istic times for recon- 
nection, iFisk fc SchwadronI (|2001[ ) have determined the 
diffusion coefficients inside and outside of a CH to be 

^, respectively. 
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Figure 4. Frequency distributions of boundary displacement ve- 
locities (a), magnetic reconnection rates (b), and diffusion coefB- 
cients (c). 



4. DISCUSSION AND CONCLUSIONS 

The relocation of open magnetic field lines through 
interchange reconnection is the most widely accepted 
method to explain some of the distinctive characteristics 
of CHs - namely their rigid rotation, the lack of loops 
with heights >10 Mm inside CHs, and the CH bound- 
ary relocations through interaction with QS loops. In 
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this study CH boundary locations were determined using 
an intensity threshold i ng tec hnique described in detail 
in iKrista fc Gallagheil (|2OO90 . This method determines 
the CH boundaries based on observations rather than al- 
lowing the observer to choose a subjective threshold. A 
number of CHs were tracked for several days to study the 
changes in the CH boundaries. The CH boundary track- 
ing was carried out by temporally coupling the bound- 
ary contours, which allowed relative distances of bound- 
ary movements and hence velocities to be established. 
In our study the displacement velocities were found to 
be isotropic with unsigned values typically <2 km s^^, 
which imply reconnection rates of <3 x 10"'^. 

These boundary velocities are in good agreement with 
the isotropic nature and ~ l-2 km s~^ av erage velocity of 
random granular motions (jMarkovlll991| ). Hence, we sug- 
gest that the majority of the reconnection at CH bound- 
aries is likely to be initiated through random granular 
motions. Velocities in excess of ~2 km s~^ were occa- 
sionally observed during our analysis. Such surges in the 
boundary displacement could be explained by consecu- 
tive interchange reconnections in the case of favourable 
loop alignments, whereby a number of loops have op- 
posite magnetic polarity to the CH flux in their nearest 
footpoint to the CH. In close proximity the CH open field 
lines could continuously reconnect with the loops and ad- 
vance into the QS, before being halted by loops that face 
the CH with the same polarity in their nearest footpoint 
to the CH. 

Interch ange reconnection has been previously dis- 
cussed bv lFisk &: SchwadronI ()2001f ) as a mechanism that 
allows the transfer of magnetic flux from one point to 
another (see Equation [ij . We determine typical values 
of diffusion coefficient to be <3 x 10^^ cm^ s~^ using 
Equation [5] with every individual CH boundary displace- 
ment measurement (i.e., from every 5 degree segment 
of each CH time pairing). Although the values derived 
in this study are comparable to those determined by 
iFisk fc SchwadronI (3.5 x lO^^ cm^ s'^ inside a CH), it 
is important to note that the results are achieved by dif- 
ferent methods. Those presented here are based on CH 
boundary tracking using observations and reflect upon 
the d iffusion processes h appening at the CH boundary, 
while IFisk fc SchwadronI provided estimates of the mag- 
netic diffusion in CHs based on typical loop heights, foot- 
point distances and sizes of supergranular cells. 

Previous studies found that the potential fleld source 
surface and magnetohydrodynamic models are success- 
ful in reconstructing the evolution of the global mag- 
netic field through a sequence of steady-state solutions. 
Despite the quasi-steady and interchange reconnection 
models disagreeing about the existence of disconnected 
open fields, it is agreed that the highly dynamic short- 
term changes in the small-scale magnetic field may be 
explained by interchange reconnecti on between open and 
closed fields (|Antiochos et al.ll2007l ). This could also be 
an alternative approach to explain how open field can 
be conglomerated (CH birth) or diffused (CH decay) on 
short-timescales. The local loop-emergence rate and the 
specific orientation of QS loops near the CH boundary 
may also influence the growth and decay of small CH 
candidates. Magnetic bipoles emerging within a CH or 
small loops admitted into the CH during interchange re- 



connection could accumulate over a CH lifetime and ulti- 
mately lead to the fragmentation and diffusion of the CH. 
Bipoles and open flux regions would have to be tracked in 
time to provide observational evidence for the processes 
governing the birth and the decay of CHs. This could also 
provide information on how small CHs turn into large, 
stable and long-lived CHs, and what processes initiate 
the decay of a long-lived CHs. 
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